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Bistability of Slow and Fast Traveling Waves in Fluid Mixtures 
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The appearence of a new type of fast nonlinear traveling wave states in binary fluid convection 
with increasing Soret effect is elucidated and the parameter range of their bistability with the com¬ 
mon slower ones is evaluated numerically. The bifurcation behavior and the significantly different 
spatiotemporal properties of the different wave states - e.g. frequency, flow structure, and con¬ 
centration distribution - are determined and related to each other and to a convenient measure of 
their nonlinearity. This allows to derive a limit for the applicability of small amplitude expansions. 
Additionally an universal scaling behavior of frequencies and mixing properties is found. 

PACS: 47.20.-k, 47.10.+g, 47.20.Ky 


In binary fluid mixtures heated from below there is an 
interesting feed-back loop between the fields of concen¬ 
tration, velocity, and temperature: The buoyancy force 
that drives the convective flow is changed by concentra¬ 
tion variations. They in turn are produced via the Soret 
effect [Q,|| by temperature gradients, i.e., via thermod¬ 
iffusion and reduced by dissipative concentration diffu¬ 
sion and by mixing through the convective flow. This 
coupling chain causes a surprising richness of spatiotem¬ 
poral pattern formation || even close to the onset of 
convection. In particular there are convective structures 
@® consisting of coupled traveling waves (TWs) of ve¬ 
locity, temperature, and concentration with significantly 
different shapes 00 - Since nonlinear concentration 
advection is typically much larger than linear diffusion 
- the ratio of these two transport rates can easily be 
above 1000 - it is not surprising that these TW states 
are typically strongly nonlinear. 

In this letter we elucidate how with increasing Soret 
coupling there appear two different bistable TWs — one 
about twice as fast as the other — which both stably co¬ 
exist with the stable quiescent conductive state. The con¬ 
vective amplitude of the fast (slow) TW is small (large) 
while the amplitude of its concentration contrast is large 
(small). The fast stable TWs have so far remained un¬ 
noticed in experiments j| |l^| and numerical simulations 
|^,^5|. They develop with increasing Soret coupling via 
a saddle node bifurcation out of a dent in the subcrit- 
ically bifurcating unstable TW branch. They are most 
easily accessible via a two-loop hysteresis at larger but 
experimentally realizable Soret effects. Furthermore, we 
discovered an universal scaling behavior of TW frequen¬ 
cies and mixing properties. And we found that the ratio 
of flow and phase velocity is a relevant parameter and 
a convenient measure of their nonlinearity that allows 
to determine where an amplitude expansion around the 
onset breaks down. 

To investigate roll like convection structures in a hor¬ 
izontal layer of, say, ethanol-water with Lewis number 
L = 0.01 and Prandtl number a = 10 we determine the 
convective solutions that bifurcate with a lateral period¬ 


icity length A = 2 out of the quiescent heat conduct¬ 
ing state 161. A finite difference method 0 as well as 
a many mode Galerkin scheme is used to solve the ap¬ 
propriate field equations 00 in a vertical cross section 
through the rolls perpendicular to their axes. Horizontal 
boundaries at top and bottom, 2 = ±1/2, are no slip, 
perfectly heat conducting, and impermeable. Our con¬ 
trol parameter r = R/R° c measuring the thermal driving 
is the Rayleigh number R reduced by the critical one 
R° c = 1707.762 for onset of convection in a pure fluid. 
We are interested here in negative separation ratios ip 
||. Then the Soret coupling between deviations T of 
temperature and C of concentration from their means 
tends to increase (decrease) the ethanol concentration in 
the cold (warm) fluid regions. 

In Fig. [I] we show how an increasing Soret coupling 
strength changes the bifurcation diagrams of (a) maxi¬ 
mal vertical flow intensity w^ lQX , (b) TW frequency w, 
(c) mixing number 1 — M, and (d) convective contribu¬ 
tion N — 1 to the Nusselt number. The order parameter 
1 -m = i-vmiicLj measures the convective 
mixing. It is defined such as to vanish in the conductive 
state, C con d = —ipz, and it approaches ±1 for convection 
with strong mixing properties where the spatially aver¬ 
aged square of the concentration variation < C 2 > be¬ 
comes small. The overall subcritical bifurcation topology 
is caused by the interplay of two adverse effects: (i) The 
Soret coupling to the degrees of freedom of the concen¬ 
tration field stabilizes the quiescent basic state [|I9| , since 
the Soret induced concentration distribution reduces the 
buoyancy force that drives convection. This shifts with 
increasing \ip\ the TW bifurcation threshold r osc (arrows 
in Fig. |l|) upwards along the r axis, (ii) With increasing 
convection mixing advectively reduces the Soret induced 
concentration gradients and with it the influence of the 
Soret effect on the buoyancy so that the convection be¬ 
havior of the mixture approaches that of the pure fluid 
dashed curves labelled ip = 0 in Fig. Qa, d. 

We found that the spatiotemporal properties of TW 
states substantially change when the convective flow 
Wmax becomes larger than the phase velocity v = 
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FIG. 1. Evolution of TW-bifurcation diagrams with Soret 
coupling strength ip: (a) squared maximal vertical flow 
(b) frequency u), (c) mixing number 1 — M, and (d) convective 
contribution to the Nusselt number N — 1 vs reduced Rayleigh 
number r. Stable (unstable) TW states are marked by filled 
(open) symbols. Four of them on the vertical dotted line at 
r = 2.24, ip = —0.6 are identified for later discussion by differ¬ 
ent triangles. Arrows mark Hopf thresholds r osc for onset of 
TW convection. The ip = 0 pure fluid limit is included in (a) 
and (d) by the dashed line. Full lines through the data points 
of (a) represent the fit discussed in the text. Only states in 
the shaded region of (a) are weakly nonlinear (cf. text). 

at r osc with a large Hopf frequency u>h and w max = 0 
the frequency (Fig. 0b) monotonically decreases along 
the TW solution branch while the convection amplitude 
Wmax grows. States with \ = w max /v < 1 (shaded re¬ 
gion in Fig. 0a) are weakly nonlinear, while those with 
larger \ are strongly nonlinear. In the former there are 
only open streamlines in the frame comoving with the 
TW (Fig. |a) and the concentration wave profile is basi¬ 
cally harmonic (Fig. 0b V and ▼). On the other hand, 
for x > 1 there are also regions of closed streamlines 
(Fig- 1 c). With growing x their size and with it the an- 
harmonicity of the concentration wave (Fig. [|b, A and A) 
increases, since within regions of closed streamlines con¬ 
centration is diffusively mixed to alternatingly high and 
low plateau levels. When moving along a TW bifurcation 
branch one observes with growing \ first an increase of 


the amplitude of a harmonic concentration wave up to 
X — 1 and then a decrease in amplitude combined with 
its wave profile becoming more 



FIG. 2. (b) Wave profiles of concentration deviation C 
(scaled by A Ta//3) at the vertical position marked by arrows 
in (a) und (c). Therein streamlines in the frame comoving 
with the TW are shown. Symbols identify states in Fig. |l] on 
the dotted vertical line. 

and more trapezoidal (Fig. 0b)- 

The structural changes that occur with growing x can 
be observed experimentally in topview and even better in 
sideview shadowgraph intensity distributions ]2G| , [To| , |T2[ . 
The sideviews in Fig. 0 of four TWs (marked by trian¬ 
gles in Fig. [l]) existing at the same r ss 2.24 show in 
(a) the smooth distribution of a weakly nonlinear linear 
state (V) and in (b)-(d) strongly nonlinear ones. In the 
latter the Soret induced ethanol rich (poor) boundary 
layers near the top (bottom) plate increasingly discharge 
via concentration jets into the growing roll like regions 
of closed streamlines. Therein, the concentration is then 
homogenized on alternating high and low levels. 

With increasing Soret coupling strength \ip\ the TW 
solution branches in Fig. [l| become more and more con¬ 
torted. Already at 'tp ~ —0.001 the lower unstable 
branches of w^ ax and N have two inflection points with 
a dent in between that can clearly be seen in Fig. [I] at 
ip = —0.25. The significance of such a dent on the lower 
unstable branch of solutions being unaccessible to exper¬ 
iments [l3|] as well as to earlier numerical simulations 
00 has not been appreciated appropriately although 
it can be seen in the paper of Bensimon et al. [ plj in 
the limit of small Soret coupling and a —> oo. With 
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our new finite difference method and our many mode 
Galerkin scheme we can now trace out the complete so¬ 
lution branch and determine how its bifurcation topology 
changes with ip. 



FIG. 3. Sideview shadowgraph intensity distributions ob¬ 
tained as described in 1 1 1(| I a| from the numerically determined 
TW states that are marked by the respective triangles in 

Fig. §■ 


We found that the above mentioned dent develops at 
xp ~ —0.4 into a forwards bending arc when the first 
inflection point (* in Fig. [|) emits two new saddles (■ 
and • in Fig. |]). So in the shaded region of the con¬ 
trol parameter plane of Fig. [| three stable states occur: 
For xp < —0.4 two different stable TWs coexist — one 
fast, the other one slow — that compete with each other 
and with the stable quiescent basic state. Furthermore, a 
new scenario appears for xp < —0.61. There, the Rayleigh 
number of the saddle with small N, w max and large to 
drops below the saddle with large N, w max and small 
u) so that a two-step hysteresis opens up between basic 
state, slow TWs, and fast TWs [ [ 22 | . 

The first inflection point (* in Fig. [|) of the TW bifur¬ 
cation branch coincides closely with the boundary \ = 1 
between weakly and strongly nonlinear states — the 
change in curvature of the bifurcation branch reflects the 


change in the TW structure. Even for xp as small as about 
— 0.001 amplitude equations of quintic order might possi¬ 
bly reproduce realistically only those TW states that lie 
between r osc and the first inflection point. The fact that 
X — 1 is an upper limit to amplitude expansions around 



¥ 

FIG. 4. xp-e phase diagram of TW states at subcritical 
driving. In the shaded parameter regime two stable TWs 
exist. Bifurcation diagrams of IV vs t in the insets explain 
the meaning of the symbols: •, ■, and ♦ represent saddle 
node states with small, medium, and large convection ampli¬ 
tude (large, medium, and small frequency), respectively. The 
first (second) inflection point is denoted by * (+). 

the Hopf bifurcation threshold follows also directly from 
the behavior of the multivalued functional dependence of 
the flow intensity A^ = Wmax on r - Its singlevalued in¬ 
verse r{X) — to see it just turn Fig. [l|a by 90° — can be 
represented extremely well by the product 

r(X)=r 0 (X)Q(X) (1) 

cf. the full lines in Fig. |Ta. Here ro(A) is the ip = 0 bor¬ 
der curve for the pure fluid (dashed line in Fig. jl]a) that 
starts at ro(A = 0) = 1 and grows for large A' slightly 
sublinearly. The rational function 

p/ _ «o + oiAT + a 2 X 2 + (I 3 A 3 
^ l + b 1 X + b 2 X 2 + b 3 X 3 

ensures with a 0 = r osc that r( X = 0 ) = r osc and with fi¬ 
nite 03/63 that r — r 0 does not diverge for large X. With 
the third order polynomials in Q four different X can 
have the same r. In that sense is a minimal represen¬ 
tation of the bifurcation diagrams of Fig. [lja. Now, the 
radius of convergence X c of an amplitude expansion of Q 
in a power series in A is given by the absolute value of 
that pole of Q closest to A^ = 0 in the complex A'-plane. 
Our fits for Q(X) show that this pole lies on the negative 
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real X-axis with X c close to v 2 so that \ = w max /v ~ 1 
is indeed an upper limit for the applicability of amplitude 
expansions. 

A recent TW model |23| into which enters X as a rel¬ 
evant parameter predicts that the frequencies are uni¬ 
versally determined by the ’’distance” r(A') — ro(X) of 
the 



1 -co/co H 

FIG. 5. Universal, ip independent scaling relations be¬ 
tween: (a) Frequency and ’’distance” r — ro of TW states from 
ip = 0 convection (cf. text), (b) degree of convective mixing 
1 — M and TW frequency. Here u>h is the Hopf frequency 
at onset. Shown are all TW states of Fig. |lj and some that 
have been suppressed there for clarity (stable TWs - filled 
symbols, unstable TWs - open symbols). 


TW states from the pure fluid convection state (if) = 0). 
For any TW with a particular A' = w^ax this distance 
parallel to the r-axis can be read off in Fig. |a. We 
confirm this prediction in Fig. |]a. It shows that the 
scaled frequencies 1 — lo/ujh are indeed linearly related 
to \Jr — ro /ujh hr an universal, ^-independent manner 
that holds for stable as well as for unstable TWs. Fur¬ 
thermore, according to Fig. ^b the convective mixing 
1 — M grows linearly with the reduced frequency devia¬ 
tion 1 — ui /ujh from the Hopf frequency and also 1 — M 
is roughly linearly related to y/r — tq/loh- 

To summarize, we have found bistability between slow 
and fast TWs and we have elucidated how the latter ap¬ 
pear when the Soret coupling becomes sufficiently strong. 
Test runs with smaller Lewis and/or larger Prandtl num¬ 
bers lead to similar bifurcation behavior. To check our 
results experimentally — preferably in annular containers 

and also to investigate effects related to spatially con¬ 
fined states, we suggest to use mixtures with ip — —0.6 
as in Refs. jf^. 
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